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Abstract 

The moments of the photon spectrum in the inclusive B — > X s 7 decay can be 
calculated to order a s accuracy at present, without knowing the a s corrections 
to the effective Hamiltonian. We discuss the standard model predictions, 
and how the moments of the photon spectrum are related to certain matrix 
elements of the heavy quark effective theory. The sensitivity of these moments 
to new physics is small, and they provide a model-independent determination 
of the b quark pole mass (at order a s ), or equivalently, the matrix elements 
A and Ai of the effective theory. 



*Work supported in part by the U.S. Dept. of Energy under Grant no. DE-FG03-92-ER 40701. 
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I. INTRODUCTION 



The inclusive B — > X s 7 decay has received a lot of attention in recent years |T]-§|, 
primarily due to its sensitivity to physics beyond the standard model (SM) @J3,0]. As any 
flavor changing neutral current process, it can only arise at one-loop level in the SM, and 
therefore possible new physics can yield comparable contributions. However, the recent 
CLEO measurement pj excludes large deviations from the SM. 

Since the h quark is heavy compared to the QCD scale, the inclusive B — > X s 7 decay 
rate can be calculated in a systematic QCD-based expansion M. The decay rate computed 
in the m& — > 00 limit coincides with the free quark decay result. Corrections can then be 
included in an expansion in powers of 1/mj, and a s (m&). 

At present, the theoretical prediction for the decay rate suffers from large uncertain- 
ties, as the result is only known in the leading logarithmic approximation. To refine the 
theoretical prediction and thus increase the sensitivity to new physics, a next-to-leading 
order calculation is needed. That is a very demanding task, as it requires the evaluation of 
many two-loop and even the infinite parts of three-loop diagrams. In the absence of such 
a calculation, and since the recent CLEO result shows no evidence for new physics con- 
tributing significantly to the B — > A s 7 decay, we investigate what the presently available 
measurement could teach us. 

We point out that the moments of the photon spectrum can be obtained to order a s 
accuracy by a relatively simple calculation. We evaluate these corrections to the first few 
moments of the photon spectrum. Since B — > X s 7 is a two-body decay at the quark level 
(at leading order in a s ), the photon spectrum is monochromatic in the spectator model. 
Therefore, the moments are also sensitive to the nonperturbative corrections in the heavy 
quark expansion. They provide a model-independent determination of the b quark pole 
mass, i.e., the matrix elements A and Ai of the heavy quark effective theory (HQET). 
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II. THE EFFECTIVE HAMILTONIAN 



The B — > X s 7 decay in the standard model is mediated by penguin diagrams. The QCD 
corrections to this process form a power series in the parameter a s hx(Myy/ml), that is too 
large to provide a reliable expansion. Therefore, it is convenient to integrate out the virtual 
top quark and W boson effects (and possible new physics) at the W scale, and sum up the 
large logarithms using the operator product expansion and the renormalization group. We 
work with the operator basis and effective Hamiltonian of Ref. || 

AC 8 

H« = --^ V* V tb Ciin) O.di) , (2.1) 

i=l 

where 

0\ = (c L/3 7^ b La ) (s La 7 M c Lf3 ) , 2 = (c La 7^ b La ) (s Lj3 7 M c L/3 ) , ^ 
7 = m b s La b Ra , 8 = ^ m b s La G% b R p . 

We listed here only the operators whose Wilson coefficients are of order unity; C3 — Cq are 
about an order of magnitude smaller, since they arise only due to operator mixing. As the 
matrix elements of all operators except O7 contain an overall factor of a s (once we use the 
"effective" Wilson coefficients || defined below), we shall neglect 3 — O e . 

The order a s corrections to the the photon spectrum come from three sources [j5|: 

(i) Corrections to Ci(M\y) coming from the matching of the SM matrix elements onto the 
effective Hamiltonian: the order a s corrections to C 7 (M W ) and C 8 (M W ) are not yet 
known. 

(ii) Corrections to the running of the Wilson coefficients C\ between the W and the m b 
scale: the next-to-leading order mixing of the dimension six with the dimension five 
operators resulting from three- loop diagrams is not known at present. 

(iii) The a s corrections to the matrix elements of the operators in the effective Hamiltonian 
at the low scale: these corrections also involve unknown two-loop diagrams. 

Certain parts of the next-to-leading log anomalous dimension matrix (ii) have been calcu- 
lated |Q. The unknown contributions are expected to be significant, as these terms should 



reduce the large /i-dependence of the leading log result. The last class of corrections (hi) 
has been considered in Ref. 0. However, important two-loop diagrams involving 2 , that 
contribute to the spectrum near the maximal photon energy, were neglected. Moreover, 
in the absence of the next-to-leading order result for CV(/V), any calculation of the spec- 
trum inevitably contains renormalization scheme dependence at order a s , and is therefore 
ill-defmed.g 

Since none of these three ingredients of the full order a s calculation are known, it may be 
surprising that the moments can be calculated to this accuracy. We prove this in the next 
section. Thus, for our purposes it is sufficient to summarize the results for the coefficients 
of the effective Hamiltonian in the leading logarithmic approximation, for which we adopt 
the scheme independent definitions of 0. In the standard model C 2 (Mw) = 1, and 

„ , , r , 3x 3 - 2x 2 . -8x 3 - 5x 2 + Ix 
° 7{Mw) = 4(^T)T lnX + 24(x - 1)3 ' 

Q -7 ' 2 ™ 3 _i_ , 2 _i o 

where 2 = mf/M^. At a low scale /1, these coefficients become 

CM = \ (V /23 + v- 12/23 ) C 2 (M W ) , 

C 7 cff (/i) = V 6/23 ^r(M^) + 5 _ ^ c&{Mw) + £ ^ ^ ? 

5 

Cf(p) = V 14 ' 23 C 8 (M W ) + C 2 (M W ) J29iV h , (2-4) 

i=l 

where 77 = a a (Mw)/a a (/i), and the numerical values of the Ttj's, g^'s, Oj's, and b^s can be 
found in ||. The scale is usually chosen to be /1 = m^, and one estimates the uncertainties 
related to the unknown higher order terms by varying typically between rrib/2 and Im^. 
However, our results will be scale independent to order a s , so this is not going to be a large 
uncertainty. 



*We thank Mark Wise for pointing this out to us. 
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III. MOMENTS 



In this section we show that moments of the photon spectrum can be calculated to 
order a s , although none of the three previously described ingredients of the full order a s 
computation of the decay rate are known. Since experimentally one needs to make a lower 
cut on the photon energy, we define the moments of the photon spectrum as 



M n (E ) 



E i dr 

Eo ^7 

E r x dr 



7 



(3.1) 



Here £™ ax 



E ^7 

[M| — (Mk + M n ) 2 ]/2Ms is the maximal possible photon energy. To illustrate 
our argument, we denote schematically the contribution of a given operator Oj to the n-th 
moment of the spectrum by (Oi) n . Then we can rewrite the moments M n as 



Cf (0 7 ) n + Cf (0 8 )n + C 2 (0 2 )r 



left 



Cf(O 7 ) + Cf (O 8 ) + C 2 (O 2 ) 



(Or)n + T^fT (Os)n + 7^ W, 



C4 



cf 



/~<en c 



(O 7 )o + T^f (O 8 ) + 



(3.2) 



For i ^ 7, the contributions (Oj)„ are of order a s . Therefore, to determine the order a s 
corrections to M n , it is consistent to take into account these contributions, and the a s 
correction to the matrix element of the operator Oj as well. But it is sufficient to know the 
Wilson coefficients Cj only to the presently available (leading log) accuracy. 

The other important observation is that the two-loop part of the order a s contributions 
to matrix elements of the operators at the low energy scale do not contribute to the moments 
of the photon spectrum. The reason is that these are virtual corrections that yield finite 
delta-function contributions at the maximal photon energy, and therefore they contribute 



equally to the numerator and the denominator of eq. ( |3.1| ). Thus, these contributions also 
cancel to order a s in the moments M n . 

We mentioned earlier that, in the absence of the next-to- leading order result for the 
coefficient C 7 , any calculation of the photon spectrum (in particular that of [0) is renor- 
malization scheme dependent at order a s . However, this scheme dependence affects again 



only the finite part of the delta-function at the maximal photon energy and therefore drops 
out from the moments. 

These arguments prove that the moments M n , as calculated below to order a s , are 
renormalization scheme and scale independent. 



IV. ORDER a s QCD CORRECTIONS 

As we mentioned in the introduction, the heavy quark expansion proves that the free 
quark decay model result is the leading term in a systematic expansion in powers of I /nib, 
in which the first nonperturbative corrections arise at order l/m\. In this section we discuss 
the order a s (rrib) corrections to the moments of the photon spectrum in the free quark decay 
model. Following ||, we keep the strange quark mass finite to regularize collinear divergen- 
cies. We use the MS subtraction scheme, Feynman gauge, and dimensional regularization 
for infrared and ultraviolet divergencies, and phase-space integrals as well. We verified the 
calculation of Ref. |9[], and also computed terms that were neglected in that calculation. The 
quantity that is simple to calculate in perturbation theory is 

Sm n (x ) = Ao pl , (4.1) 

dT v 



FQDM 



dx 



dx 



where r pQDM denotes the decay rate in free quark decay model, and we introduced the 
dimensionless parameters 

x = - , r = -| . (4.2) 

(1 — r) mi, nit 

The variable x corresponds to E 1 j _E™ ax in the free quark decay model. The definition ( |4.1| ) 

makes it apparent that the functions 5m n (xo) are proportional to a s and that they are not 

affected by corrections proportional to 5(1 — x) at this order. In the last section we shall 

discuss how to relate 5m n to the experimentally measurable moments M n . 

The conclusion of Ref. |J is that near the maximal photon energy (x = 1) only the 

operators O2 and O7 are important. Although we found that the interference of O7 and Og 



is also peaked near x = 1, this term turns out to be numerically small fllO|| . In this letter we 
only include the dominant contributions that were already calculated in M, while we shall 



present the full order a s calculation of the moments M n elsewhere |10[ . 

The contributions involving the operator O2 are regular functions of the photon energy, 
and their explicit form can be found in |5|. The only singular (and numerically the most 
significant) contribution to the the free quark decay rate r FQDM at order ot s comes from the 
operator O7 alone. In the r — > limit this contribution reads (we only explicitly present in 
this letter the corrections in the r — > limit, but we included the r-dependent terms in our 
numerical results): 

dr 77 



dx 



~m h {n)- 


2 { 







+ 



a s C F 



4tt 



1-^(5 + 1^-21*^ 
4tt V 3 /i ; 



7 + x - 2x 2 - 2(1 + x) ln(l - x) 
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5(1 -x) 

ln(l -x) 



(4.3) 



+ 4- 



x 



X 



where Cp = 4/3 in 577(3), and 



G 2 F \V tb VZ\ 2 aCf(fif 

32 7T 4 



ml 



(4.4) 



By m b we mean the b quark pole mass at order a s , as discussed below. The [/(x)] + distri- 
bution corresponding to a function f(x) acts on a test function g(x) as 

[f(x)]+g(x)dx= [ f(x)[g(x)-g(l)}dx. (4.5) 



'0 ^0 

By including the /z-dependence in eq. (|4.3|) , we can explicitly verify that T 77 is /i-independent 
to order a s . This cancellation by itself does not reduce significantly the /z-dependence of the 
theoretical prediction for the total decay rate, as the dominant part of C 7 efr (/z) arises from 
the mixing of O2 with 7 at order a s , but still at leading log. The a s correction in eq. ( |4.3|) 
modifies the prediction for the total decay rate by about 15%, while it affects the first two 
moments of the photon spectrum, Mi and M 2 , by less than 3% and 5%, respectively. Our 
numerical results for 5m\ and 5m 2 , including the contribution of O2, are shown in Table I. 

When evaluating these corrections, it has to be kept in mind that the perturbative 
expansion becomes singular in the photon endpoint region, and a resummation of the per- 
turbative corrections may be required. As we are interested in the moments M n for small n, 



and (a s CV/27r) ln 2 (l — x), the exponent of the Sudakov factor that suppresses the endpoint 
spectrum becomes of order unity only around x ~ 0.99, our calculation is consistent 
without taking these effects into account. 



V. NONPERTURBATIVE CORRECTIONS 

We include the nonperturbative corrections to the free quark decay to order l/m 2 and 
leading order in a s . At that order only corrections to the matrix element of the operator 
O7 contribute. We present here the resulting corrections to the photon spectrum Jl2[ only 
in the r — > limit again: 



The dimensionful constants Ai and A2 parameterize the matrix elements of the kinetic and 
chromomagnetic operators, respectively, which appear in the Lagrangian of the HQET at 
order 1/toq: 

X 1 = ^(M(v)\h v (iD) 2 h v \M(v)) , 

A 2 = -j- (M{v) I f h v <j af3 K \M(v)) , (5.2) 

where dp = 3 and dy = —1 for pseudoscalar and vector mesons, respectively. M(v ) denotes 
the meson state and h v denotes the quark field of the effective theory with velocity v. The 
numerical value of A2 can be extracted from the mass splitting between the vector and 
pseudoscalar mesons, A 2 = (m|, - m|)/4 ~ 0.12 GeV 2 , while there is no similarly simple 
way to determine Ai from experiments (see, e.g., |13|])- Without including the order a s 
corrections discussed in the previous section, one obtains to order 1/m 2 in the heavy quark 
expansion 

m b / Ai + 3A 2 \ M B - A 



Ai + 3A2 r//1 A Ai 



2m 2 

<*?>-W a = ~ ■ (5-3) 



In general, the central moments {{E 1 — {E 1 )) n ) for n > 2 are proportional to Ai(mb/2) n ~ 2 . 
Therefore, they are particularly useful for measuring Ai. This is not unexpected, since Ai 
is the measure of the Fermi motion of the b quark that is responsible for the smearing of 
the photon spectrum. A comparison of the values of Ai extracted from different central 
moments can be used to estimate the systematic errors. 



The parameter A in eq. (|5.3|) describes the mass difference between a heavy meson and 
the heavy quark that it contains, and it is one of the parameters that set the scale of the 
1/m expansion jl4| . It is related to the B meson mass and the b quark pole mass according 



to 



M B = mb + A-±±lh. (5.4) 

IVTLh 



The quantity A suffers from renormalon ambiguities ||15|| . However, at any finite order in a s 
it is consistent to extract A (or the pole mass nib) from certain experiment (s), and use the 
resulting numerical values to evaluate theoretical predictions accurate to the same order in 
perturbation theory for other processes 

The series of the nonperturbative corrections to the moments of the photon spectrum, 
M n (Eo), is under control if the invariant mass of the final hadronic state corresponding to 
the lower cut Eq is above the resonance region. Eq. ( |5.1|) is only related to the experimentally 
measured spectrum once the theoretical expression is smeared over typical hadronic scales. 
This smearing is provided by taking the moments M n (E ), if n is not too large, and Eq is 
sufficiently far from E™ x . Using the relation 



E 7 , (5.5) 



2M B 

we see that the present experimental signal region of E 1 > 2.2 GeV corresponds to 
Mx s ~ 2.2 GeV. Even below this scale the widths of the X s resonances are typically larger 
than their mass differences, and there are no resonances above 2.5 GeV [|Tj]]. It is still 
important for the reliability of our analysis to try to lower the experimental cut on the photon 
energy. For example, E = 2 GeV or Eq = 1.8 GeV would correspond to Mx s ^ 2.6 GeV or 
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Mx s ~ 3 GeV, respectively. Varying E provides a check on the systematic uncertainties: 
the extracted values of A and Ai should be unaffected by the variations of E , once the 
corresponding hadronic invariant mass is sufficiently above the resonance region. 



VI. SUMMARY AND CONCLUSIONS 

We can summarize our discussion by writing the theoretical prediction for the first two 
moments of the photon spectrum as 



M B - A 



/2E \ / 2 A 2 A 3 



M 1 (E ) = ^— - l + 5m 1 (^)+Ola 2 s ,a s ^ 1 ,^) , (6.1) 

Z L \ mi, ' ^ m b m b - 

We used A to denote some QCD scale of order Aq CD or A. To the order these relations are 
accurate, it is consistent to replace mj by Mb — A everywhere in eq. ( |6.1| ). The numerical 
results for Smi(xo) and 5rri2(xo) are listed in Table I for three different values of xq. 

The significance of these relations is that they provide a reliable means of determining 



A (that is, A at order a s |1CJ) and Ai, or equivalently, measure the b quark pole mass at 
order a s . Especially the first relation in eq. ( |6.1|) is remarkable, since it is independent of 
Ai (and A2), and very sensitive to A, with small theoretical uncertainties. The sensitivity of 
the second relation to Ai is reduced because of the factor 1/12. 

We would like to emphasize that the left-hand side of the relations in eq. (|6.1| ) are 
measurable at CLEO; in fact, the central values can be extracted from Ref. [0]. As we do 
not know the cross-correlations of the errors on the data points, we are not in a position to 
quote numerical values for the experimental uncertainties. From the central values of the 
data, solving the first equation in ( |6.1| ), we find (with large uncertainties) 

A~ 450MeV, m b ~ 4.83 GeV. (6.2) 

We do not quote even a central value for Ai, as the present experimental data do not 
constrain it to any reasonable accuracy. By varying all input parameters {Cj{fi) and C^/i) 
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corresponding to the range m&/2 < \i < 2mb, a s corresponding to 0.11 < a s {Mz) < 0.13, 
m s between 100 MeV and 500 MeV, m t , and m c ) we find that the theoretical uncertainty of 
this measurement of A will be as small as ±30 MeV, while that of Ai about ±0.15GeV 2 . 

In view of our earlier discussion, it is important to try to expand the experimental signal 
region. On the one hand, the systematic uncertainties inherent in our analysis (related to 
how well duality holds) can be estimated by varying the lower cut on the photon energy, as 
discussed at the end of Section V. On the other hand, expanding the signal region would 
diminish the sensitivity of the results as to whether the Sudakov logarithms at the endpoint 
are resummed or not. 

The sensitivity of the moments of the spectrum to new physics is limited by how much 
operators other than O-j affect M n . We found that O2 does not contribute to 5mi and 8iri2 
by more than 10% in the SM. Given that the experimental constraint on the total decay 
rate from CLEO excludes large deviations from the SM, we conclude that the moments 
are largely insensitive to new physics. Even if physics beyond the SM contributes to the 
B — > X s 7 decay, the proposed determination of A and Ai is likely to remain unaffected. 

We conclude that the moments of the photon spectrum in the inclusive B — > X s 7 
decay will provide reliable measurements of fundamental parameters of QCD, which in 
turn will refine theoretical predictions for other observables in heavy quark decays. The 
complete order a s calculation of the moments and a more detailed analysis of the theoretical 
uncertainties will be presented in a forthcoming paper [|K] . 
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TABLES 





x = 0.91 


x = 0.83 


x = 0.75 


Smi(x ) 


r = 4 • 10" 4 


-0.014 


-0.020 


-0.025 




r = 1 • 10~ 2 


-0.012 


-0.017 


-0.020 


5m 2 (xo) 


r = 4 • 1(T 4 


-0.028 


-0.040 


-0.046 




r = 1 • 1(T 2 


-0.023 


-0.032 


-0.038 



TABLE I. Central values of 5mi(xo) and 5rri2(xo) for two different values of m s . r = 4 • 10 4 
corresponds to m s = 100 MeV, while r = 1 • 10~ 2 corresponds to a constituent quark mass 
m s = 500 MeV. For m b ~ 4.8 GeV, x = 0.91, x = 0.83 and x = 0.75 correspond to E = 2.2 GeV, 
E = 2.0 GeV and = 1.8 GeV, respectively. 
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